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Abstract. Numerous theoretical studies using various equation of state models have shown that quark 
matter may exist at the extreme densities in the cores of high-mass neutron stars. It has also been shown 
that a phase transition from hadronic matter to quark matter would result in an extended mixed phase 
region that would segregate phases by net charge to minimize the total energy of the phase, leading to the 
formation of a crystalline lattice. The existence of quark matter in the core of a neutron star may have 
significant consequences for its thermal evolution, which for thousands of years is facilitated primarily 
by neutrino emission. In this work we investigate the effect a crystalline quark-hadron mixed phase can 
have on the neutrino emissivity from the core. To this end we calculate the equation of state using the 
relativistic mean-field approximation to model hadronic matter and a nonlocal extension of the three- 
flavor Nambu-Jona-Lasinio model for quark matter. Next we determine the extent of the quark-hadron 
mixed phase and its crystalline structure using the Glendenning construction, allowing for the formation 
of spherical blob, rod, and slab rare phase geometries. Finally we calculate the neutrino emissivity due 
to electron-lattice interactions utilizing the formalism developed for the analogous process in neutron star 
crusts. We find that the contribution to the neutrino emissivity due to the presence of a crystalline quark- 
hadron mixed phase is substantial compared to other mechanisms at fairly low temperatures (< 10 9 K) 
and quark fractions (< 30%), and that contributions due to lattice vibrations are insignificant compared 
to static-lattice contributions. 

PACS. 21.65.Qr quark matter - 26.60.Kp equations of state - 97.60.Jd neutron stars 


1 Introduction 

Seconds after the formation of a neutron star the mean 
free path of neutrinos grows beyond the star’s radius, and 
neutrinos from the core escape easily cooling the new star 
rapidly. Neutrino emission continues to be the dominant 
energy loss mechanism of a neutron star for thousands of 
years until a temperature of about 10 7 K is reached pQ. It 
has been shown that the presence of quark matter in the 
core of a neutron star can have a significant impact on the 
neutrino emissivity, and suggested that this impact could 
have an observable effect on the star’s thermal evolution 
m- In this work we investigate the effect that a mixed 
phase containing quark matter may have on the neutrino 
emissivity of a neutron star. 

If electric charge neutrality in a neutron star is to 
be treated globally as proposed by Glendenning j4], then 
the first order phase transition from hadronic matter to 
quark matter in the core will result in a mixed phase in 


which both phases of matter coexist. To minimize the to¬ 
tal isospin asymmetry energy the two phases will segre¬ 
gate themselves, resulting in positively charged regions of 
hadronic matter and negatively charged regions of quark 
matter, with the rare phase occupying sites on a Coulomb 
lattice. Further, the competition between the Coulomb 
and surface energy densities will cause the matter to ar¬ 
range itself into energy minimizing geometric configura¬ 
tions [2|- 

The presence of the Coulomb lattice and the nature 
of the geometric configurations of matter in the quark- 
hadron mixed phase may have a significant effect on the 
neutrino emissivity from the core. More specifically, the 
presence of electrons in the mixed phase will lead to an ad¬ 
ditional neutrino emissivity mechanism due to interactions 
with the lattice. This process is analogous to neutrino- 
pair bremsstrahlung of electrons in the neutron star crust, 
where ions exist on a lattice immersed in an electron gas, 
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and for which there exists a large body of work (see, for ex¬ 
ample [51[6ll7ll8l[9l l 10 imn ). The situation is somewhat more 
complicated in the quark-hadron mixed phase, but the op¬ 
erative interaction is still the Coulomb interaction. Thus, 
to estimate the neutrino-pair bremsstrahlung of electrons 
in the quark-hadron mixed phase we rely heavily on this 
body of work (particularly ini'). 

Neutrino emissivity due to the interaction of electrons 
with a crystalline quark-hadron mixed phase has been pre¬ 
viously studied in this manner by Na et al. m- In the 
present work we replace the MIT Bag Model used by m 
to describe quark matter with a three-flavor nonlocal vari¬ 
ant of the Nambu-Jona-Lasinio model. Next, we extend 
the range of possible geometric structures in the mixed 
phase beyond spherical blobs to include rods and slabs, 
and calculate the associated static lattice contributions 
to the neutrino emissivity. Phonon contributions to the 
emissivity for rod and slab geometries are not considered, 
though a comparison of the phonon and static lattice con¬ 
tributions for spherical blobs is given and indicates that 
phonon contributions may not be significant. Finally, the 
extent of the conversion to quark matter in the core is de¬ 
termined for the chosen parameterizations, and this allows 
for a targeted comparison between emissivity contribu¬ 
tions from standard neutrino emission mechanisms (mod¬ 
ified Urea, nucleon-nucleon and quark-quark (NN+QQ) 
bremsstrahlung) and contributions from electron-lattice 
interactions. In this work the minimal cooling paradigm 
is assumed, as the mechanism under investigation is not 
expected to compete with the direct Urea process, but 
may serve to enhance the cooling of neutron stars in its 
absence. 

This paper is structured as follows. In Sec. [2] we dis¬ 
cuss the zero temperature equation of state of a neu¬ 
tron star containing hadronic matter, quark matter, and 
a quark-hadron mixed phase. The crystalline structure 
of the mixed phase is described in Sec. [3] The neutrino 
emissivity due to interactions between electrons and the 
crystalline lattice in the quark-hadron mixed phase is de¬ 
scribed in Sec. 21 Our results, including neutron star prop¬ 
erties and neutrino emissivity calculations, are presented 
in Sec. [3] Finally, we present our conclusions in Sec. [6] 


2 Equation of State 


2.1 Neutron Star Crust 


The neutron star outer and inner nuclear crust exists at 
densities between 10 4 g cm -3 < ecrust ^ 10 14 g cm -3 
m- Matter in the inner crust consists mostly of nuclei in a 
Coulomb lattice that is immersed in a gas of electrons and, 
above neutron drip (> 4 x 10 11 g cm -3 ), free neutrons. 
In this work we use a combination of the Baym-Pethick- 
Sutherland and Baym-Bethe-Pethick equations of state for 
the nuclear crust 1X811791 . Compared to the neutron star 
core the crust has less affect on the neutron star properties 
that are to be studied in this work. 


2.2 Confined Hadronic Phase 


The hadronic phase of neutron star matter exists at den¬ 
sities above that of the crust and is populated by baryons 
(B = {n,p,A,£,£}) and leptons (A = {e~,p~}). To 
model the hadronic phase we use the relativistic mean- 
field approximation (RMF), in which the interactions be¬ 
tween baryons are described by the exchange of scalar (er), 
vector (w), and isovector (p) mesons |20| . The mean-field 
Lagrangian is given by [T51[TTllT51l211l22ll25] 


C=^2i/j b 


B 


7m ( id>L - ^ 9 p bt ' 


- (m n - g a Bcr) 


0 b + -{d^ad^a - m 2 a a 2 ) 


- -b, 


3 1 _ / _ _\4 


- 


3 -am n ( g„<j ) - -c a (g a cr) 

1 

4 C 

1 

2 to p^m ' P'' ~ 4 ■ P r 

Y 0A (*7m^ - rn\) 0A ■ 


( 1 ) 


The er and u> mesons are responsible for nuclear binding 
while the p meson is required to obtain the correct value 
for the empirical symmetry energy. In contrast to er and ui 
mesons, which are isoscalars, the p meson is an isovector 
field that manifests itself in the occurrence of the Pauli 
matrix r (= (t 1 ,? - 2 ,r 3 )) in Eq. (fTjl. The cubic and quar- 
tic er terms in Eq. © are necessary (at the relativistic 
mean-field level) to obtain the empirical incompressibility 
of nuclear matter mm- The field tensors and p pv 
are defined as uj pv = 9^-9^ and p pv = d p p v -d v p p . 

The meson-baryon coupling constants (g a N, 9uN, gpN , 
b a , c a ) of the Lagrangian are set so that the properties 
of nuclear matter at saturation density are reproduced 
for the appropriate parameterization (Table [1]). In this 
work we employ the GM1 and NL3 parameterizations as 
in Ref. [23]. To fix the meson-hyperon coupling constants 
g m y we follow the method presented in Ref. [33]. The 
scalar meson-hyperon coupling constants g a y are fit to the 
following hypernuclear potentials at saturation density: 
U ( a N) = -28 MeV, U { S N) = +30 MeV, and ui N) = -18 
MeV. The vector meson-hyperon coupling constants g u y 
are fixed in SU(3) flavor symmetry by the mixing angle 
0 V and coupling ratio z taken from the Nijmegen extended- 
soft-core (ESC08) model [35] . The isovector meson-hyperon 
coupling constants g p y are given by the usual relations, 
g p A = 0 and g p s = 2 g p ~ = 2g P N- 

The field equations for the baryon fields follow from 
Eq. © as follows [21 IT3inXHI3] . 

(n^M ~ m s) 0s = -gaBCTlpB + g^Bl^^pIpB 

+9pBl M T ■ p^lpB ■ (2) 
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Table 1. Properties of nuclear matter at saturation density for 
the NL3 and GM1 parameterizations. Properties include the 
nuclear saturation density po, energy per nucleon E/N, com¬ 
pression modulus K, effective nucleon mass m * N , and asymme¬ 
try energy a sy • 


the states’ chemical potential. The baryonic and leptonic 
number densities ( ps , p\) are both given by 

p i = (2J i + l)-^. (12) 


Nuclear saturation properties 

GM1 

NL3 

CcE 

1 

O 

Cl 

0.153 

0.148 

E/N (MeV) 

-16.3 

-16.3 

K (MeV) 

300 

272 

m*/m n 

0.78 

0.60 

a sy (MeV) 

32.5 

37.4 


The meson fields in © are solutions of the following field 
equations mmmm, 

(<9^ + ml)a = — m N b N g a N (ScrAfCr) 2 

B 

-C N g a N (gcrNCr) , (3) 

+ m^u„ = ^ g^B^Blv^B , (4) 

B 

+ ™ 2 pPv = Yl dpB^BT'yulpB ■ (5) 

B 


The free parameters of the system are the meson mean- 
fields (<r, lo, p), and the neutron and electron fermi mo¬ 
menta (k n , fc e ). Finally, the energy density and pressure 
of the hadronic phase are given by mmm 


£h = ^ bm n (g a a ) 3 + ^c(g a a) 4 + 7 } m l<j 2 

+ 


E 


2Jb + 1 
27T 2 


J \Jk 2 + m* B (a) k 2 dk 


+ Ei/ \/k 2 +m 2 x k 2 dk , 


(13) 


PH = - Tbm n (g G af - ^ c(g a a) 4 - ^m 2 cr 2 
+ 0 + o m pP03 


In the mean-field limit, the meson field equations ([3]) through 
([5]) are given by pi flMmiTh] 


2 \ ' 2 Jb 

m a a = 2 ^ g*B 

B 


r^B 




27T 2 


xfk 2 




(o') 


k 2 dk 


( 6 ) 



- bm n g a (g (7 a) 2 - cg a (g a a) 3 , 

E guiB 

ia pb ' (7) 

B ^ 

P03 = ^ B pB , (8) 

B m P 

where the effective baryon mass m* B (<j) = ms — g<jBa- 
To determine the equation of state we solve a nonlin¬ 
ear system consisting of the meson mean-field equations 
and the charge conservation conditions (baryonic, electric) 
given by jJMSQSj 

Pb - E P B = 0 > ( 9 ) 

B 

E PBVB +]T P\Q\ = 0 , (10) 

B X 


2.3 Deconfined Quark Phase 

If the dense interior of a neutron star contains deconfined 
quark matter, it will be made of up («), down (d), and 
strange (s) quarks in chemical equilibrium with a small 
number of electrons and muons. To model the quark phase 
we use a nonlocal extension of the Nambu-Jona-Lasinio 
model (n3NJL) as described in Ref. [25]. The effective 
action of this model is given by 

Se = J d 4 x{ip(x)(ifl-rh)i/j(x) + i Gs[ (ip(x)X a ip(x)) 2 

+ (VS(x)f 7 5 A a V’(x)) 2 ] + H [det[ip(x)(l + j 5 )^(x)] 
+det[V>(x)(l - 75)V’( a 0] ] - Gy[(^(a;)7 / "A a '0(a:)) 2 
+{ip(x)i'y li 'y 5 \a*lj(x)) 2 } j, (15) 


where pb is the total baryonic density and qs and q x are 
the electric charges of baryons and leptons, respectively. 
Particles in the hadronic phase are subject to the chemical 
equilibrium condition, 

Pi = hp N - qip e , (11) 

where pi is the chemical potential and bi is the baryon 
number of particle i. New baryon or lepton states are pop¬ 
ulated when the right side of equation CD is greater than 


where / denotes quark flavor ( u,d,s ), ^ is a chiral U(3) 
vector that includes the light quark fields, ip = (u, d , s) T , 
rh = diag(m u , mj, mn s ) is the current quark mass matrix, 
A a with a = 1,8 denote the generators of SU(3), and 
Ao = ^/2/3 113 x 3 - The coupling constants Gs and H , the 
strange quark mass m s , and the three-momentum ultravi¬ 
olet cutoff parameter A, are all model parameters. Their 
values are taken from Ref. [26] , i.e., m u = rrid = 5.5 
MeV, m s = 140.7 MeV, A = 602.3 MeV, G S A 2 = 3.67 
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and HA 5 = — 12.36. The vector coupling constant Gy is 
treated as a free parameter. 

For the mean-field approximation, the thermodynamic 
potential associated with Se of Eq. m is given by 


and the charge conservation equations, 

E Pf - 3ph = 0 , (24) 

f—u,d,s 


flNL = -? E . 

f=u,d,s 


E 

f=u,d,s 1 

E 

- f~ U :d,S 

of 2 


dp 0 / dp 



dpp [(pf -E f )0(p f -m f )\ 


~ T f S f 


G S c2 
~2~ bf 




S d S s 


f 


^' 4Gy’ 

f—U,d,S 


(16) 

where cry, wy, and Sf are the quark scalar, vector, and 
auxiliary mean fields, respectively. Moreover, we have Ef = 

^jp 2 + Wy, Wy = (p 0 + i-Pf) 2 + P 2 i and My(wy) = my + 

(7fg(uj y) are the momentum dependent quark masses. The 
quantity g{u) 2 ) = exp(— uj 2 / A 2 ) is the form factor which 
introduces nonlocality into the quark interactions m- 
The auxiliary mean fields are given by 


Sf = -48 / dp 0 


dp , 2x 


M f(u 2 f) 

UJ 1 + Mf(ldf) 


(17) 


E p f q f + E pxqx = 0 • ( 25 ) 

f=u,d,s \=e~ ,n~ 

Finally, the pressure and energy density are given by 

p Q = n o - 12 NL , (26) 

e Q = -PQ + E PfPf + E PxPx > ( 27 ) 

f=u,d,s 

where l?o is the grand thermodynamic potential fI NL cal¬ 
culated for pf = uJf = 0. 


2.4 Quark-Hadron Mixed Phase 

When the pressure in the hadronic phase grows to a level 
equal to that of the quark phase at the same baryonic 
density a first order phase transition from hadronic mat¬ 
ter to quark matter may begin. Since a theory that can 
treat both the hadronic and quark phases simultaneously 
is currently unavailable, we construct the mixed phase by 
blending RMF and n3NJL. Each phase is solved for sep¬ 
arately, and then the two are blended together under the 
Gibbs condition, pn = Pq ■ The pressure (jpm) and energy 
density (cm) in the mixed phase are given by |3]|4] 

Pm = ^(ph+Pq), (28) 


Due to the inclusion of the vector interaction the quark 
chemical potentials are shifted as follows, 


Pf = Pf ~ g(™ 2 f)uf , 


(18) 


Wy = (po + ipf) 2 +p 2 ■ (19) 

The scalar and vector mean fields are obtained by mini¬ 
mizing the grand thermodynamic potential, 


<9l? NL 

da f 


= 0, 


<912 NL 

duff 


= 0. 


( 20 ) 


The quark number densities are given by 


<912 NL 

Pf dp f 


( 21 ) 


To determine the equation of state one must solve a 
nonlinear system of equations for the fields efy and <Dy, 
and the neutron and electron chemical potentials p n and 
p e . This system of equations consists of the mean field 
equations, 

di + GsSi + —HSjSk = 0 , ( 22 ) 

with cyclic permutations over the quark flavors, 


ujf - 2 Gy d = 0 , (23) 


and 

cm = (1 - x)ch + XtQ , (29) 

where x = Vq/V Total is the quark fraction of the mixed 
phase. Other properties such as the particle number den¬ 
sities can be handled in a similar fashion. 


3 Crystalline Structure of the Quark-Hadron 
Mixed Phase 

A mixed phase of hadronic and quark matter will arrange 
itself so as to minimize the total energy of the phase. 
Under the condition of global charge neutrality this is 
the same as minimizing the contributions to the total en¬ 
ergy due to phase segregation, which includes the surface 
and Coulomb energy contributions. Expressions for the 
Coulomb (ec) and surface (eg) energy densities can be 
written as |3] 

e c = 2ne 2 [q H (x) ~ Qq{x)} 2 r 2 xf D {x) , (30) 

e s = Dxa{x)/r, (31) 

where qn (qcf) is the hadronic (quark) phase charge den¬ 
sity, and r is the radius of the rare phase structure. The 
quantities x and fo{x) in Eq. (l30l) are defined as 

x = min(x, 1 — x) 


(32) 
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Fig. 1. Equations of state of this work. Shown are the hadronic 
and mixed phases. Square markers indicate the beginning of 
the quark-hadron mixed phase, and dot markers indicate the 
location of the maximum mass. 


and 


Id(x) 


1 

D + 2 


^2 (2 


(33) 


where D is the dimensionality of the lattice. The quantity 
er(x) in Eq. (IdTll denotes the surface tension. 

The phase rearrangement process will result in the for¬ 
mation of geometrical structures of the rare phase dis¬ 
tributed in a crystalline lattice that is immersed in the 
dominant phase (Figure [2]). The rare phase structures are 
approximated for convenience as spherical blobs, rods, and 
slabs [3]. The spherical blobs occupy sites in a three di¬ 
mensional (D = 3) body centered cubic (BCC) lattice, 
the rods in a two dimensional (D = 2) triangular lattice, 
and the slabs in a simple one dimensional (D = 1) lattice 
mi- At x = 0-5 both hadronic and quark matter exist as 
slabs in the same proportion, and at x > 0-5 the hadronic 
phase becomes the rare phase with its geometry evolving 
in reverse order (from slabs to rods to blobs). 



Quark Quark Quark Hadronic Hadronic Hadronic 

Blobs Rods Slabs Slabs Rods Blobs 


0 «- V Q /V -► 1 


constructing a single theory that can accurately describe 
both hadronic matter and quark matter. Therefore, we 
employ an approximation proposed by Gibbs where the 
surface tension is taken to be proportional to the difference 
in the energy densities of the interacting phases [3] , 

cr(x) = r]L [e Q (x) - e H {x)] > (34) 

where L is proportional to the surface thickness which 
should be on the order of the strong interaction (1 fm), and 
rj is a proportionality constant. In this work we maintain 
the energy density proportionality but set the parameter 
rj so that the surface tension falls below 50 MeV fm -2 , a 
value consistent with those suggested for cr(x) in recent 
literature l35ll36113711381 . 


3.2 Rare Phase Structure Size, Charge, and Number 
Density 


The size of the rare phase structures is given by the ra¬ 
dius (r) and is determined by minimizing the sum of the 
Coulomb and surface energies, d ^ ec d ^ es ' > , and solving for r 

0' 


Dcr(x) 


4tt e 2 f D {x) [qH (x) - Qq(x)Y 


(35) 


The primitive cell of the lattice is taken to be the Wigner- 
Seitz cell, though it is simplified to have the same geom¬ 
etry as the rare phase structure. The Wigner-Seitz cell 
radius R is set so that the cell is charge neutral. 

The density of electrons in the mixed phase is taken to 
be uniform throughout. Charge densities in both the rare 
and dominant phases are also taken to be uniform, an 
approximation supported by a recent study by Yasutake 
et al. [35] . The uniformity of charge in the rare phase also 
justifies the use of the nuclear form factor ( F(q )) presented 
in Section U The total charge number per unit volume 
(\Z\ /Vr a re) of the rare phase structures is given in Figure 

SI 

The number density of rare phase blobs will be impor¬ 
tant for calculating the phonon contribution to the emis¬ 
sivity. Since there is one rare phase blob per Wigner-Seitz 
cell, the number density of rare phase blobs (rib) is simply 
the reciprocal of the Wigner-Seitz cell volume, 


n b = (4 ttR 3 /3)~ 1 . (36) 


4 Neutrino Emissivity in the Quark-Hadron 
Mixed Phase 


Fig. 2. (Color online) Spherical blob, rod, and slab rare phase 
structures. Vq/V Total denotes the quark fraction of the mixed 
phase. 


3.1 Surface Tension of the Quark-Hadron Interface 

Direct determination of the surface tension of the quark- 
hadron interface is problematic because of difficulties in 


Modeling the complex interactions of electrons with a back¬ 
ground of neutrons, protons, hyperons, muons, and quarks 
is an exceptionally complicated problem. However, to make 
a determination of the neutrino emissivity that is due 
to electron-lattice interactions in the quark-hadron mixed 
phase we need only consider the Coulomb interaction be¬ 
tween them. This simplifies the problem greatly, as a sig¬ 
nificant body of work exists for the analogous process of 
electron-ion scattering that takes place in the crusts of 
neutron stars. 
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Fig. 3. (Color online) Rare phase structure radius r and Wigner-Seitz cell radius R plotted against quark fraction for the given 
parameterizations and vector coupling constant values. The change in background color refers to a change in structure as shown 
in the right panel. Discontinuities in the radii are also associated with changes in the crystalline structure. 



Fig. 4. (Color online) Rare phase structure charge number 
per unit volume plotted against quark fraction for the given 
parameterizations and vector coupling constant values. Slight 
discontinuity at \ = 0.5 is due to a higher number of posi¬ 
tively charged baryons (electrons make up difference in nega¬ 
tive charge). 


4.1 Electron-Lattice Interaction 


To determine the state of the lattice in the quark-hadron 
mixed phase we use the dimensionless ion coupling param¬ 
eter given by [l] 


Z 2 e 2 




RkbT 


(37) 


Below T me it = 175 the lattice behaves as a Coulomb liq¬ 
uid, and above as a Coulomb crystal pQ. It was shown in 
Na et al. [12] that the emissivity due to electron-blob in¬ 
teractions in the mixed phase was insignificant compared 
to other contributions at temperatures above T > 10 10 



Fig. 5. Temperature below which Umklapp processes are 
frozen out (Tumkiapp) as a function of quark fraction, and con¬ 
tributions to the neutrino emissivity due to electron-phonon 
interactions become negligible. 


K. Therefore, in this work we consider temperatures in 
the range 10 7 K < T < 10 10 K. At these temperatures the 
value of the ion coupling parameter is generally well above 
Tmeitj and so the lattice in the quark-hadron mixed phase 
is taken to be a Coulomb crystal. 

To account for the fact that the elasticity of scattering 
events is temperature dependent we need to compute the 
Debye-Waller factor, which is known for spherical blobs 
only and requires the plasma frequency and temperature 
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given by 


CjJ v — 


/ 4:irZ 2 e 2 rib 
m b 


(38) 


= m 

where m b is the mass of a spherical blob El- The Debye- 
Waller factor is then given by 


W(q) 


otq 

8fej 


(1.399 ( 


-9.lip 


12.972 tp) spherical blobs, 
rods and slabs, 


(40) 

where q = |g| is a phonon or scattering wave vector, 
a = AH 2 k^/(kBT p m b ), and t p = T/T p [TTII341 . In order 
to smooth out the charge distribution over the radial ex¬ 
tent of the rare phase structure we adopt the nuclear form 
factor given in m, 




3 

W) 


[sin(gi?) — qRcos(qR)] . 


(41) 


Screening of the Coulomb potential by electrons is taken 
into account by the static dielectric factor e(q, 0 ) = e(q), 
given in Ref. [28]. However, the charge number of the rare 
phase structures is high and the electron number density 
is low, so setting this factor to unity has no noticeable 
effect on the calculated neutrino emissivity. Finally, the 
effective interaction is given by m 


processes in which a phonon is created (or absorbed) by an 
electron that is simultaneously Bragg reflected, resulting 
in a scattering vector q that lies outside the first Brillouin 
zone, q 0 > ( 67 r 2 ^) 1 / 3 |171I30| . where n b is given by Eq. 
(1561) . 

The contribution to the neutrino emissivity due to 
phonons is contained in L pb/ and given by Eq. (21) in Ref. 

m, 


L p h = 


r 1 s eS (qmg )\ 2 
y 0 V y\e(qM 2 



(45) 


where y = q/(2k e ), and the lower integration limit y p ex¬ 
cludes momentum transfers inside the first Brillouin zone. 
The structure factor S e g is given by Eqs. (24) and (25) in 

Ref. fITj). 


S eS (q) = 189 - 

7 r 


-2 w 


da 


1 - 40£ 2 + 80£ 4 
(1 + 4 £ 2 ) 5 cosh 2 ( 7 r£) 

x (e*«> - l) , 

(46) 

m = _ 008 M _\ (47) 

2 m b \ oj s sinh (Hu s /2k bT) / ’ 

where a = tksT/h and (...) denotes averaging over phonon 
frequencies and modes, 


</-(*)> 



(48) 


V(q) 


ATrepzFjq) w(q) 
q 2 e(q) 6 


(42) 


It is assumed that there are three phonon modes s, two 
linear transverse and one longitudinal. The frequencies of 
the transverse modes are given by = a^fc, where i = 


4.2 Neutrino Emissivity 

General expressions for the neutrino emissivity due to 
electron-lattice interactions were derived by Haensel et al. 
m for spherical blobs and by Pethick et al. EH for rods 
and slabs, 

Qbiobs « 3.35 x 10 - 67 n b T 6 Z 2 L McV s _1 fm -3 , (43) 

Q rods,slabs » 3.00 x 10 " 88 k e T 8 J MeVs " 1 fm " 3 , (44) 

where L and J are dimensionless quantities that scale the 
emissivities. Both L and J contain a contribution due to 
the static lattice (Bragg scattering), but we consider the 
additional contribution from lattice vibrations (phonons) 
for spherical blobs, so L = L s \ + L pb . 

4.3 Phonon Contribution to Neutrino Emissivity 

The expressions for determining the neutrino emissivity 
due to interactions between electrons and lattice vibra¬ 
tions (phonons) in a Coulomb crystal, with proper treat¬ 
ment of multi-phonon processes, were obtained by Baiko et 
al. 133] and simplified by Kaminker et al. . The phonon 
contribution to the emissivity is primarily due to Umklapp 



Fig. 6. (Color online) Mass-radius relationship of neutron stars 
for the given parameterizations and vector coupling constant 
values listed in Table [T] Square markers indicate the beginning 
of the quark-hadron mixed phase, and dot markers indicate 
the location of the maximum mass. The green shaded region 
indicates the mass constraint set by PSR J3048+0432. 

1, 2, cii = 0.58273, and 02 = 0.32296. The frequency of the 
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longitudinal mode coi is determined by Kohn’s sum rule, 

m . 

Umklapp processes proceed as long as the temperature 
lumkiapp ^ T p Z 1 / 3 e 2 /(He), below which electrons can no 
longer be treated in the free electron approximation [30 . 
This limits the phonon contribution to the neutrino emis¬ 
sivity to only a very small range in temperature for a 
crystalline quark-hadron mixed phase (see Figure [5]), and 
renders it negligible compared to the static lattice contri¬ 
bution as will be shown in the next section. 


4.4 Static Lattice Contribution to Neutrino Emissivity 


Pethick and Thorsson m found that with proper han¬ 
dling of electron band-structure effects the static lattice 
contribution to the neutrino emissivity in a Coulomb crys¬ 
tal was significantly reduced compared to calculations per¬ 
formed in the free electron approximation. Kaminker et 
al. m presented simplified expressions for calculating the 
static lattice contribution (L s \) using the formalism devel¬ 
oped in Ref. 110] . The dimensionless quantities L s \ and J 
that scale the neutrino emissivities for spherical blobs and 
rods/slabs, respectively, are given by 


F s i 


and 


i v (i-^)IfxaT 

12 Z^ 0 y\ \e(I<)\ 2 


I(y K ,tv)e- 2W ^ 

(49) 



(50) 


where K = \K\ is a scattering vector and restricted to 
linear combinations of reciprocal lattice vectors, yx = 
K/(2k e ), t v = k B T/ [\V(K)\(1 - y 2 K )\, and I(y K ,t v ) is 
given by Eq. (39) in Ref. [H]. The sum over K in Eqs. (HU) 
and m terminates when K > 2 k e , prohibiting scattering 
vectors that lie outside the electron Fermi surface. 


5 Results 

The neutron star equation of state has been calculated us¬ 
ing the relativistic mean field approximation to describe 
the hadronic phase and the three-flavor nonlocal Nambu- 
Jona-Lasinio model for the quark phase, with the Gibbs 
condition governing the combination of the two in the 
mixed phase (Figure [l]). Using the equation of state we 
solve the Tolman-Oppenheimer-Volkoff equation mm 
and find the mass-radius relationships given in Figure 0 
The maximum masses of the neutron stars obtained for 
the given parameter sets are able to account for the re¬ 
cently discovered high mass pulsars PSR J3048+0432 and 
PSR J1614-2230 [4211431144] . excluding GM1 with no vec¬ 
tor coupling. It is evident from Figure [G] that increasing 
the vector coupling constant increases the maximum mass 
for the particular parameterization. 

Figures [7] and [5] show the relative particle densities for 
the NL3 and GM1 parameterizations and three different 



PjPo 


Fig. 7. (Color online) Relative particle number densities as 
a function of baryonic density (in units of nuclear saturation 
density) for the NL3 parameterization and given values of the 
vector coupling constant. Grey shaded region indicates den¬ 
sities beyond the maximum mass neutron star for the given 
parameterization. 



pJpq 


Fig. 8. (Color online) Same as Figure [7] but for the GM1 
parameterization. 


values of the quark vector coupling constant. Hyperoniza- 
tion does not occur at all in the NL3 parameterization, as 
it is preceded by the low density onset of the quark-hadron 
phase transition at 2 — 3 times nuclear density. The same 
is true of the GM1 parameterization except in the case 
that Gy = 0.10 Gs- Here the onset of the quark-hadron 
phase transition occurs at a much higher density due to 
the presence of the S~ and A hyperons which soften the 
equation of state considerably, an effect that can be seen 
in the right panel of Figure [1] The low density onset of 
the quark-hadron phase transition is due in part to the 
choice of meson-hyperon coupling constants, which have 
been shown to postpone the onset of hyperonization, stiff¬ 
ening the low density equation of state [51]. Figure |5] shows 
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Table 2. Properties of the maximum mass neutron star for the given parameterizations. 


Gv/Gs 


GM1 



NL3 


0 

0.05 

0.10 

0 

0.05 

0.10 

M/Mq 

1.89 

2.05 

2.20 

2.04 

2.24 

2.43 

Xmax 

0.32 

0.31 

0.15 

0.31 

0.30 

0.32 

Pb [l/fm 3 ] 

0.75 

0.76 

0.80 

0.61 

0.61 

0.62 

e [MeV/fm 3 ] 

851 

883 

969 

687 

696 

726 


that neutron stars within about 0.1 —0.2 M@ of their max¬ 
imum mass contain a quark-hadron mixed phase in their 
core, with most possessing a maximum quark fraction of 
around 30% (see Table [2]). 

Figure [10] shows the neutrino emissivity that is due to 
the crystalline structure of the quark-hadron mixed phase 
for all parameterizations and temperatures between 10' — 
10 10 K, as well as the modified Urea and bremsstrahlung 
(NN+QQ) emissivities for comparison. Electron-phonon 



Fig. 9. (Color online) Comparison of the static lattice and 
phonon contributions to the neutrino emissivity for the NL3 
parameter set at T > Tumkiapp- The comparison is performed 
for only the spherical blob geometry, as the phonon contribu¬ 
tion is not determined for other geometries. 

interactions contribute to the neutrino emissivity when 
the mixed phase consists of spherical blobs (x < 0-21 and 
X > 0.79) and only at T > T Umk i app (Figure 0. Figure El 
shows that the static-lattice contribution to the emissivity 
dominates the phonon contribution rendering it negligi¬ 
ble, particularly at quark fractions relevant to the neutron 
stars of this work (x < 0.3). Equations (TT51) and (I5U1) in¬ 
dicate that the static-lattice contribution is calculated as 
a sum over scattering vectors Nk that satisfy K < 2 k e . 
At the onset of the mixed phase the electron Fermi mo¬ 
mentum k e is at a maximum, which is particularly large 
in magnitude due to the lack of hyperons that would typ¬ 
ically aid in the charge neutralization process. However, 
as the quark-hadron phase transition proceeds the neg¬ 


atively charged down and strange quarks take over the 
process of charge neutralization, resulting in a rapidly de¬ 
creasing electron number density ( k e = (37 t 2 n e ) 3 ). This 
and the exponentially decreasing size of the Wigner Seitz 
cell in the spherical blob phase (x < 0.21) lead to the steep 
decline in Nk (Figure [TT]), which accounts for the rapid 
decrease of the neutrino emissivity in the mixed phase. 

The geometrical structure of the quark-hadron mixed 
phase terminates with the rod phase at \ ~ 0.3 for nearly 
all the chosen parameterizations. Up to this point the neu¬ 
trino emissivity due to the structure of the mixed phase 
is either larger or comparable to the modified Urea and 
bremsstrahlung (NN+QQ) emissivities for T < 10 8 K. 
The emissivities in the NL3 and GM1 parameterizations 
are comparable, though the effect of the mixed phase struc¬ 
ture appears more substantial for NL3 due to lower modi¬ 
fied Urea and bremsstrahlung (NN+QQ) emissivities. The 
emissivity at very low quark fraction (x < 0.05) may be 
overestimated due to the finite blob radius at x = 0 that 
results from the fact that fd( 0) = 0.4. Finally, beyond 
the rod-slab structure transition at x ~ 0.35 the electron- 
lattice contribution to the overall neutrino emissivity is 
negligible, though this is beyond the extent of the mixed 
phase of the neutron stars in this work. 


6 Summary and Conclusions 

Exploring the properties of compressed baryonic matter, 
or, more generally, strongly interacting matter at high 
densities and/or temperatures, has become a forefront area 
of modern physics m ■ Experimentally, such matter is be¬ 
ing created in relativistic particle colliders such as the Rel¬ 
ativistic Heavy Ion Collider RHIC at Brookhaven and the 
Large Hadron Collider (LHC) at Cern, and great advances 
in our understanding of such matter are expected from the 
next generation of collision experiments at FAIR (Facil¬ 
ity for Antiproton and Ion Research at GSI) and NICA 
(Nucloton-based Ion Collider fAcility at JINR) [151150] . 

Complementary to these experiments, astrophysics pro¬ 
vides a natural laboratory in which to explore the physics 
of compressed baryonic matter too (see EUE1M! and ref¬ 
erences therein). The Hubble Space Telescope and X-ray 
satellites such as Chandra and XMM-Newton in particular 
have proven especially valuable. New astrophysical instru¬ 
ments such as the Five hundred meter Aperture Spheri¬ 
cal Telescope (FAST), the square kilometer Array (skA), 
Fermi Gamma-ray Space Telescope (formerly GLAST), 
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Fig. 10. (Color online) Comparison of the neutrino emissivity due to electron-lattice interactions in the quark-hadron mixed 
phase to modified Urea and bremsstrahlung (NN+QQ) processes in the mixed phase for different parameterizations. Shading 
represents geometric structures as shown in Figure [3] 



Fig. 11. (Color online) The number of scattering vectors that satisfy the condition K < 2 k e as a function of the quark fraction 
in the mixed phase for the equations of state of this work. 


Astrosat, ATHENA (Advanced Telescope for High EN- 
ergy Astrophysics), and the Neutron Star Interior Com¬ 
position Explorer (NICER) promise the discovery of tens 
of thousands of new neutron stars. Of particular interest 
will be the proposed NICER mission, which is dedicated to 
the study of the extraordinary gravitational, electromag¬ 
netic, and nuclear-physics environments embodied by neu¬ 
tron stars. NICER will explore the exotic states of matter 
in the core regions of neutron stars, confronting nuclear 
theory with unique observational constraints. 

With that in mind, we focus in this paper on quark 
deconfinement in the cores of neutron stars. The neutron 
star equation of state for cold catalyzed matter (T < 1 
MeV) has been determined using the relativistic mean 


field (RMF) approximation to model the hadronic phase 
and the nonlocal three-flavor Nambu-Jona-Lasinio model 
(n3NJL) for the quark phase. The mass-radius results in¬ 
dicate that a neutron star containing quark matter in the 
core can account for the high mass of the recently discov¬ 
ered pulsars PSR J3048+0432 and PSR J1614-2230, and 
that a maximum mass neutron star can be expected to 
contain approximately 30% quark matter at the center. If 
the surface tension between hadronic and quark matter is 
low as suggested in the recent literature, a phase transition 
that results in a mixed phase will occur in the core of a 
neutron star. The relaxed condition of global charge neu¬ 
trality will lead to charge segregation in the mixed phase 
resulting in the formation of a crystalline lattice of quark 
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matter immersed in a hadronic matter background. Ex¬ 
panding on Na et al. (12] . we considered the presence of 
two additional geometrical structures in the mixed phase 
in addition to spherical blobs: rods, and slabs (Figure [2]) ■ 
Using the formalism developed for analogous neutrino- 
pair bremsstrahlung processes in the neutron star crust we 
have estimated the neutrino emissivity due to electron- 
lattice interactions in the quark-hadron mixed phase. The 
emissivity is highly dependent on the electron number 
density, which has been shown to decrease considerably in 
the presence of negatively charged hyperons and quarks 
(Figures [7] and [5]) ■ However, we have shown that at tem¬ 
peratures between 10 7 K and 10 9 K and quark fractions 
less than around 30% the neutrino emissivity due to electron- 
lattice interactions is significant when compared to the 
standard baryon and quark modified Urea and bremsstrahlung' 
(NN+QQ) processes fFigure HU1). Further, we have also 
shown that the emissivity due to electron-phonon inter¬ 
actions is insignificant compared to contributions from 
Bragg diffraction at temperatures above which Umklapp 
processes are frozen out (Figure [9]). 

Before we can determine the effect the presence of 
quark matter and the crystalline structure of the quark- 
hadron mixed phase has on the thermal evolution of a 
neutron star the following steps need to be taken. First, 
RMF should be replaced with a model for hadronic matter 
that softens the equation of state and produces results for 
neutron star radii that are more compatible with obser¬ 
vations and recent statistical studies (see for example S3 
146]). To this end, we are currently working on combining 
an RMF model that accounts for density dependence in 
the values of the meson-baryon coupling constants with 
the three-flavor n3NJL model (see for example [551f571f551 
sramM. Next, the thermal conductivity and spe¬ 
cific heat should be calculated for the quark-hadron mixed 
phase as outlined in Na et al. m using the updated equa¬ 
tion of state for quark matter (n3NJL) and accounting 
for additional rare phase geometries (rods, slabs). Finally, 
these results would be incorporated into a neutron star 
cooling simulation capable of properly accounting for the 
complexity of the crystalline quark-hadron mixed phase. 
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